Faraday waves in quasi-one-dimensional superfluid Fermi-Bose mixtures 
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Generation of Faraday waves in superfluid Fermi-Bose mixtures in elongated traps is investigated. 
The generation of waves is achieved by periodically changing a parameter of the system in time. Two 
types of modulations of parameters are considered, first a variation of the fermion-bosons scattering 
length, and secondly the boson-boson scattering length. We predict the properties of the generated 
Faraday patterns and study the parameter regions where they can be excited. 
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I. INTRODUCTION 

Faraday waves (FW) are spatially periodic patterns 
that can be generated in a system with a periodic vari- 
ation in time of the system parameters. First Faraday 
waves was observed by Faraday for a vessel with a liq- 
uid oscillating in the vertical direction [1 . Such type 
of structures can exist in nonlinear optical systems pl- 
[6] where variations along the longitudinal direction of 
the Kerr nonlinearity can be achieved by the periodic 
variation of the effective cross-sectional area of the non- 
linear optical fiber. Recently parametric resonances in 
modulational instability of electromagnetic waves in pho- 
tonic crystal fibers with a periodically varying diameter 
has been experimentally observed in [7]. Other exam- 
ples are the patterns in the Bose-Einstein condensates 
(BEC) with an atomic scattering length or a radial con- 
finement (a transverse frequency of the trap) periodically 
varying in time. For BECs the cases of single [8 -10 and 
two components condensates [11] [12] , as well as dipolar 
condensates [131 IHl? have been investigated. In BEC 
FW was first predicted in 2D with a periodically varying 
transverse frequency of the trap [8 and later observed in 
an experiment with a repulsive interacting BEC, loaded 
into an elongated trap [15 . The existence of FW also in 
elongated fermionic clouds was discussed in the work [16] 
and Faraday patterns in a superfluid Fermi gas has been 
investigated recently in [TTj . 

Periodic modulation of the coefficients of nonlinearity 
in the relevant mean-field equations can be achieved by 
variation of the atomic scattering length by the Feshbach 
resonance technics [T8H2T] or by time modulation of the 
transverse frequency of the trap [HI. In the former case 
it is necessary to vary the external magnetic field in time 
near the resonant value. The presence of a deep optical 
lattice has been shown to suppress the Faraday pattern 
generation [22] . 

The purpose of this work is to investigate the mech- 
anism of Faraday wave generation in elongated super- 
fluid Fermi-Bose (FB) mixtures. Such FB mixtures have 
many interesting properties in comparison with the pure 



bosonic case [23^ and Faraday waves can be an useful 
tool to measure the nonlinear properties and in partic- 
ular instabilities in these systems. Two types of atomic 
scattering lengths are relevant in this system. First the 
fermion-bosons scattering length describing the scatter- 
ing between fermionic and bosonic atoms a/5, secondly 
the boson-boson scattering length a^. Variation in time 
of these lengths by an external magnetic field opens up for 
the possibility of generating Faraday waves in the mix- 
ture. The action of these variations are different though. 
In the latter bosonic case we parametrically excite the 
BEC subsystem, while in the first case we excite the 
bosonic and fermionic subsystems simultaneously. Thus 
we can expect different responses on temporal parametric 
perturbation with various types of pattern formation. 

Strongly repulsive interacting bosons in ID, so 
called Tonks- Girardeau gases [24], have the same long- 
wavelength dynamics as non-interacting fermions [25 . 
Hence, the results presented here can also be re- 
alised in systems with two coupled bosonic species with 
strong/weak intra-species interactions. 

II. THE MODEL 

The model of a quasi-one-dimensional superfluid FB 
mixture is described by the following system of coupled 
equations for the complex functions '0i,2(^,^) [22] 

#2,^ = -V^2,xa^ + /^7r^|V^2|'^V^2 +^2l|V^lPV^2. (1) 

In general Bose- and Fermi subsystems are described by 
the Lieb-Liniger and the Gaudin-Yang theories respec- 
tively. Here we are interested in weak Bose-Bose interac- 
tions (we consider small positive ^5) and attractive Fermi- 
Fermi interactions and the superfluid Fermi-Bose sys- 
tem is then described by the nonlinear Schrodinger-like- 
equation ([T]) [26-30 . In the BCS weak attractive cou- 
pling limit the fermionic subsystem coefficient is /t: = 1/4, 
while in the molecular unitarity limit it is z^: = 1/16 [ 23] . 
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Finally, for the bosonic Tonks- Girardeau limit [25^ with 
components 1 (2) being a weakly- (strongly-) repulsive 
bosonic species we have hi = 1. Furthermore, = 
2hai)UJ± is the ID coefficient of mean-field nonlinearity 
for bosons, where is the scattering length and ujj_ is 
the perpendicular frequency of the trap. Similarly, gi2 is 
the inter-species interaction coefficient [23 . The system 
is written in dimensionless form using the variables 



mi)UJ± 



X 



2mi)l 



Go, 



where Gj, j = 6, 12 is the coefficients of the mean- field 
nonlinearities. We also implicitly assume = m/ in ([T]), 
since such conditions can be realized approximately in 
the mixtures '^Li - ^Li and ^^K - ^^K. 

Below we will consider the case of identical interaction 
coefficients g2i = gi2- The mean- field nonlinearities will 
be varied in time as 



9l2{t) 

9b{t) 



-ai2 cos((^i2t)], 
- aj) cos{Qi}t)]. 



(2) 



Such variation can be achieved for example by using Fesh- 
bach resonance technics, namely by variation of an exter- 
nal magnetic field near a resonant value [18 . This leads 
to the temporal variation of inter- and intra-species scat- 
tering lengths and the respective mean-field coefficients. 
The system ([T]) has plane wave solutions 



where 



V^l,2 = ^l,2exp(z^l,2), Ai,2 e 



-gbAlt^gi2Alt, 



(3) 



(4) 



III. 



MODULATIONAL INSTABILITY OF 
PLANE WAVES 



equations containing Vt and Q^, we get the following sys- 
tem 

V{k,t)u^oul{t)V{k,t) = e(t)Q{k,t), 
Q{k,t)u^ujlQ{k,t) = e(t)V{k,t), 

U{k,t) = k^ J dtV, P{k,t) = J dtQ. (7) 

where 

ujl{t) = e [e + 2^5^^?] , = k\k^ + 

e{t) = 2g^2{t)k^A,A2. (8) 

Here gb{t) and gi2{t) are defined in Eqs. ([2|. In the 
following we use the notation ui = uoi(t = Ti/2Vt})) and 
£q = £(t = 11/20.12)' Hence, by solving the coupled equa- 
tions for V and Q for given /c, all the components of 
5ipi^2 can in principle be obtained by the inverse Fourier 
transform. 

We here first consider the case of a FB mixture with 
constant system parameters. Looking for solutions V 
and Q with a time dependence of the form exp(±zr^t), 
we obtain the dispersion relation of the modulations 



-4 



(9) 



Such that in the weak coupling limit £0 <^ ^2 ~ ^1 

^2 ^2 

— ^ "-^l — — 2 — ^ — 2 ' ^2 — ^ "^2 ' ^ 



UJ2 — ^\ ^2 — '^1 

while in the limit of approaching frequences UJ2 ^1 

^2 , , .2 ^ o2 V , .2 



UJo — 00- 



2 ' 



(10) 



^1 ^ ^1 ~ ^0, 1^2 ^ ^2 + ^0- 



(11) 



The stability condition (for any k) is obtained by re- 
quiring VLi^2 in (l9|) to be real, i.e. 



Let us now study the modulational instability (MI) of 
nonlinear plane waves using the linear stability analy- 
sis [5 . We will look for solutions of the form 

V^l,2 = [Ai,2+(5V^l,2(x,t)]e^^^'^('\ |(5V^1,2| « Ai,2. (5) 

Substituting these expressions into the system ([T]) and 
linearizing, we get the following system for 5^^1,2 

iSilji^t - gh{t)A\{8^l)i + 8^1)1) + S^Iji^xx 

+ ^12(^)^1^2(^2 + = 0, 

+ gi2(t)A^A2{8^i^5ri) = ^' (6) 

With difji = u -\- iv and 5tlj2 = p -\- iq^ we use the Fourier 
transform in the spatial coordinate to transform Eqs. (|6| 
into functions of the wavenumber k and the time t de- 
noted by the corresponding capital letters. Differenti- 
ation with respect to time of the remaining differential 



g['^ < ^2i:^Kg^^^Al (12) 
In opposite case we have MI in the FB mixture [[23] . 



IV. ANALYSIS OF PARAMETRICALLY 
EXCITED INSTABILITIES 

Below we will consider the linear ODE model ^ for 
different cases of periodic modulations ([2]). 

When only the inter- species interaction parameter gi2 
is modulated, we obtain a system of two coupled oscilla- 
tors with a coupling parameter varying in time 



Vtt+^lV = e{t)Q, 
Qtt+^lQ = e{t)V, 

where e{t) = eo + ai2£o cos(f2i2t) according to Eq. 



(13) 
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method [3TJ [32]. The derivations are cohected in the 
Appendix where we conclude that parametric resonances 
may occur at the two frequencies of excitations 



= 002 — ^1, ^12,+ = + ^1- 



(15) 



As the detailed analysis shows, the excitations for l^i2,- 
are stable and do not lead to FW. Inserting the expres- 
sions from Eq. (|8|) for the second case of Eq. ( 15 ), with the 



2gl^^Al and 6 



use of a simplified notation a — 
we get the resonance frequences in terms of the wavenum- 
bers of the corresponding Faraday waves 



(16) 



FIG. 1: (Color online) Solutions to the ODE model of 
Eq. ITsf . (a) Amplitudes growing exponentially in time, 
(b) Spectrum for the modulation frequency ^^i2,+ = 65.5 for 
the wavenumber k — 1 (see text) and with Qi and from 
Eq. (|9|, which agrees well with ( |1Q[ ) in this regime. Blue 
dashed curves are for V (bosons) while red thin curves are 
for Q (fermions). Initial conditions are V/Ai — Q/A2 — 10~^ 
iVt — Qt — ^) and the parameters here are: ai2 = 0.25, Ai — 



/300, A2 = V20, g^^ 



(0) 



0.01. 
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0.8 and At = 1/4. 
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FIG. 2: (Color online) Solutions to the full PDF model of 
Eq. ([1]) for driven coupling between fermions and bosons with 
the modulation frequency ^^12,+ = 65.5 {k = 1). We plot 
the oscillating parts \ — Aj) /Aj for j = 1 (bosons) in 
(a), and for j — 2 (fermions) in (b). Initial conditions are 
il^j — Aj[l + 10~^ exp(iA:x)] and boundary conditions in x are 
periodic. The parameters are the same as for Fig. [1] 

When only the intra-species parameter ^5 is modu- 
lated, we get a system of one Mathieu equation coupled 
to an oscillator equation 

Vu^ul(t)V = £oQ, 
where ojl{t) = cjf + Q^bCo'f cos{Qbt). 



(14) 



A. Driven coupling between fermions and bosons 

We consider the case of a stable FB mixture accord- 
ing to (12) such that the MI is absent [23]. Hence, any 
instability is only due to parametric resonance. To ana- 
lyze resonances of the system (13) we use the multiscale 



We invert (16), such that the wavenumbers are 



+ 6- 2^a6 + l]?2,+ 



^n2,+ - 



(17) 



The spatial wavelength L = 27r/k of the FW obtained 
from Eq. (l7\ is 



27r 



n 



12,+ ' 



a + b-2Jab + n\2,+ 



(18) 



According to the Appendix, the region of instability 
for 1^12,+ from (15) is restricted by the lines 



2 I ^0<^12 
■^1 ± 



2 I ^0<^12 /^2 



(19) 



while the exponential growth rate of the amplitudes is 
restricted by the maximal gain 



(20) 



A mathematically particularly simple case of the 
driven inter-species modulation is when uoi ^ UJ2. In- 
troducing the symmetric [^+(t) = {V -\- Q)/2)] and an- 
tisymmetric [£,-{t) = {V — Q)/2] combinations in the 
ODE (13), we can find that a set of parametric reso- 
nances at 1^12 — 21^1, 21^2 exist. This follows directly 
from the fact that we then have two uncoupled Mathieu 
equations [Sj in the variables ^+ and ^_ . Note that res- 
onances at low (21^1) and high (21^2) frequencies occur 
[see Eq. 



Numerical solutions of the ODE model (13) are pre- 
sented in Fig.[l]for parameters corresponding to the BCS 
regime. It is found that amplitudes are growing expo- 
nentially in time and the behavior agrees well with the 
theoretical predictions. Using the same parameters, we 
perform numerical simulations of the full PDF model ([T]) 
describing the Fermi-Bose mixture, see Fig. |2j Results 
shows quantitative agreement for the PDF and ODE 
models upto t ~ 10. From Fig. [3j where data for the 
spatial wavelength L of the FW versus the modulation 
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FIG. 3: (Color online) Spatial wavelengths of the Faraday 
patterns. In (a) we use parameters as for Fig. [l] with the 
left dot (k — 1) corresponding to ^^12,+ = 65.5. In (b) the 
case corresponding to Fig. [l] (right dot) is the second (blue) 
curve from the bottom. Here, as well as for the lowest (red) 
curve [left square in (a)] where ^^12,+ = 135 {k = 2), the 
condition \12\ is not fulfilled in the entire domain, hence MI 
is possible in the left part (dashed lines). The condition (12) 
for not having MI can be satisfied in the entire n range, e.g. 
by choosing A2 — 2 ■ larger, top (green) curve. Finally, 
the special case oi a — h (i.e. with g^^^ oc k) is illustrated 
with the third (black) curve, where the left dot corresponds 
to the case illustrated in Fig. |4] 
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FIG. 4: (Color online) Solutions to the ODE model of 
Eq. (14). (a) Amplitudes growing exponentially in time, 
(b) Spectrum for the modulation frequency Q^h = 63.3 {k = 1) 



and with Qi and Q2 agreeing well with (11) in this regime. 
Same initial conditions as in Fig.[l] while the parameters here 



are: ah — 0.25, Ai 
0.8 and n = 1/16. 



V5000, A2 



0.1, 



(0) 
9l2 



frequency are plotted in (a), it is seen how the wavelength 
is decreasing when the frequency is growing. While re- 
sults for the wavenumber k = 1 have been presented in 
Figs. [1] and |2j we have numerically confirmed the results 
for = 2,3,4 [squares in Fig. |3] (a)] also with full PDE 
simulations. In Fig. [3] (b) we show how L depends on 
the fermionic subsystem parameter n within the analytic 
model (p|). 



The results of Eq. ( 19 ) in practice here means that with 
the parameter values of Fig. [l] the modulation frequency 
can be in the order of 1% larger (or smaller) when excit- 



FIG. 5: (Color onhne) Solutions to the full PDE model of 
Eq. ([1]) for driven nonlinearity for bosons with the modulation 
frequency = 63.3 {k — 1). Same initial and boundary 
conditions as in Fig. |2] while the parameters are the same as 
for Fig. [4] 



ing the FW. This have been confirmed numerically with 
Eqs. ([1]) and (13), and we have also observed a weak de- 
pendence on ai2 for the maximum of the Vti2 resonance 
region. Hence, by optimizing a larger amplitude can 
be obtained [i.e. the numerical curve (A) in Fig. [6] can be 
moved further towards the line of the maximal theoretical 
gain] . 



B. Driven nonlinearity in the bosonic subsystem 

The case where only the intra-spieces (boson-boson) 
interaction parameter is modulated is described by the 
system ([l4|. Applying again the multiscale approach (see 
the Appendix), we have concluded that resonances occurs 
at VL}) = 2uJi under the additional condition that uoi ^002. 
Hence, from Eq. {Isl) with oui = 002 we now obtain 



nb = 2kVk^^a, A:^ = ^ (^yjl ^ ^l/a^ - 1^ 



(21) 



The corresponding instability region is bounded by the 
lines 



-2 2 I 



-, UJ2 CJi, 



and the maximal gain is now equal to 



4cji 



(22) 



(23) 



Numerical solutions of the ODE model (14) with ex- 
ponentially growing amplitudes are presented in Fig. [4] 
for parameters corresponding to the molecular unitarity 
limit. In the spectrum we can see peaks corresponding 
to the combination of two excited frequencies for both 
components. These predictions are again well confirmed 
by PDE simulations of the system ([T]) with driven non- 
linearity of bosons, see Fig. |5] Although the bosonic 
component is lagging behind (i.e. lower amplitude in 
Figs. [4] and [5| it is growing exponentially with the same 
rate as V^2- 



According to Eq. (22 ) the modulation frequency can be 



in the order of 3% larger or smaller (for the parameters 



5 




FIG. 6: (Color online) Exponential growth rates p of the 
slowly varying amplitudes [envelope oi V ^ Q ^ exp(pt) 
of ([7|]. Solid curves shows numerical data (dots) for the three 
cases: A (Fermi-Bose modulation, aj = ai2), B (Bose-Bose 
modulation, aj = at) and C (Super-resonance, aj = at = 
—ai2). Dashed lines shows theoretical predictions according 
to Eq. (20) for case A, and Eq. (23) for case B. In case A, 

^^12 



65.5, while in cases B (C) we use Qb — 63.3 (= Q12) 



V. GROWTH RATES OF THE FW 

In general we note, that in the regime where the lin- 
ear stability analysis based on Eqs. ^ and ([6| are 
valid we have quantitative agreement between the full 
mod el ([T| ), see Figs. [2] and |5) and the ODE models ([l3| 
and (|14|7 Fig. [l] respectively Fig. [4j We found an expo- 
nential growth in time of the oscillating amplitudes in 
this regime. For larger times the nonlinearities of the full 
model ([1]) cause a saturation of the amplitudes, while the 
results of the linear stability analyse becomes unphysical. 

In Fig. [6] we show results of the exponential growth rate 
for the three sections A, B and C above, together with 
the theoretical results for the maximal gain derived in the 
Appendix. We show results only for the wavenumber k = 
1, however, we note that the theoretical estimates Pm{k) 
are increasing with k and are asymptotically approaching 
the constant \imk ^ 00 Pm{k) = ajg^2 ^1^2/"^- 



VI. CONCLUSION 



of Fig. [4| when exciting the FW. This have been con- 
firmed numerically with Eqs. ([T]) and (14), although the 
numerical resonance region was found for slightly lower 
values of than the theoretical estimate. 

Finally we can analyse also this case in terms of the 
(anti-) symmetric combinations From Eq. (14) we 
then have the driven coupled Mathieu-like system 

^±,u + nl2^± + a.ujl cos((^5t)[^+ + ^-]/2 = 0, (24) 

and in agreement with our findings above, the literature 
here states that 21^1, 21^2 and Qi + 1^2 in this case, 
see e.g. [11133,. 



C. Driven super-resonance 



In conclusion we have illustrated the possibility of 
Faraday patterns for Fermi-Bose mixtures, i.e. with 
atomic bosons coupled to fermions, both in the fermionic 
BCS regime and in the molecular unitarity limit. In par- 
ticular we have investigated quasi-one-dimensional super- 
fluid FB mixtures with periodic variations in time of the 
Fermi-Bose or Bose-Bose interactions. We find Faraday 
patterns for both cases and study their properties de- 
pending on the parameters for modulations and of the 
system settings. Combining the two types of modula- 
tions can results in even larger amplitudes. 

We also conjecture that Faraday waves can be observed 
in an atomic BEC coupled to a Tonks- Girardeau gas. 

A natural continuation of this work is to investigate 
Faraday patterns in FB mixtures for the 2D and 3D cases. 
This problem reserves a separate investigation though, 
since the corresponding coupled nonlinear- Schrodinger- 
like equations are different. 



Comparing ( 16 ) with ( 21 ) , we note that when a = 6, we 
have 1^12,+ = ^b' Now we change a sign in Eq. ([2|, such 
that for example for ai2 = —0^5, then the two types of 
modulations initiate opposite phases for the two compo- 
nents and an increased growth rate of the Faraday waves 
are observed, as compared to the two distinct cases dis- 
cussed before. For demonstrating super-resonance nu- 
merically we used the parameters as for Fig. |4] with 
ai2 = —0.25. Also this result can be understood from 
a consideration of the system ^ for symmetric and an- 
tisymmetric combinations, ^+ and For example in 
the case when ai2 = —0^5/2 and Q12 = ^^6, we simply 
have the driven coupled Mathieu-like system 



e±,U + [^1,2 + ^2,1^ COs(l]5t)]e± 



2^5 



(25) 
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Appendix 

To understand the resonances in the coupled sys- 
tem (13), we can use the multiscale analysis [32]. Fol- 



lowing this approach, we look for solutions where 



^01 ■ 



- Soai ■ 



-ela2 - 



^02 " 



-£0^1 



-£0^2- 



(26) 
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and correspondingly for the functions V and Q 

V = Vo{t,T)^eoVi{t,T)^elV2{t,T)^..., 

Q = Qo{t,T)^eoQi{t,T)^elQ2{t,T)^..., (27) 

where T = is a slow time. Taking the terms of each 
order in eq we obtain from (13) the following system of 



equations up to the linear order in eq 
Vo,tt + ^01 "^0 = 0, Qo,tt + ^02^0 = 0, 



(28) 

-ul^Vi = -2Vo^tT -aiVo^[l^ai2Cos{nt)]Qo, 



Qi,u+^02Qi = -^Qo^tT -biQo^[l^ai2Cos{Qt)]Vo. 

The solutions of the two upper uncoupled equations 
in ( [28| ) can be written in the form 

Vo{t,T) = Ao{T)cos{uJoit)^Bo{T)sm{ujoit), 
Qo{t,T) = Co{T)cos{uJo2t)^Do{T)sm{ujo2t).{29) 

Then we require the absence of the resonant terms in 
the right hand side of the two lower equations in (28). 
With the ansatz Q = a;2 ± cji, and after averaging in 
the fast time we obtain a system of equations for the 
envelope functions Aq^Bq^Cq and I^o of (29) 





aiBo T 


2 " 


= 0, 




—2ujqiBq^t — 


aiAo + 




= 0, 




2uJo2Co,T — 


biDoT 


^Bo 
2 " 


= 0, 




— 2a;o2^o,T - 


- biCo + 


2 " 


= 0. 


(30) 



Looking for solutions of the form Aq^ Bq^Cq^ Dq ~ 
exp(pt), i.e. (for example) with Ao,t ^ P^o/^o, we find 
from ( 30 ) the characteristic equation + Mp^ -\- N = 



with coefficients 

2bluj'^^ + 2aluj'^2 T 0^12^01^02 



M = 



8^01^02 



Remember that the sign ^ in M of (31) 
002 i ^1 respectively, and note a lso that 
from -M/2±^yM^/4 - N it is seen 
can correspond to a positive real p, i.e. 
exponentially growing amplitude, while 
fti2- are stable. 

The maximal exponential growth rate 
ni2 + is found from (I3T| with - AN 



(4ai6i - 

256uj^^uj^2 

(31) 

is for fli2,± = 
> 0. Hence, 
that only l^i2,+ 
a FW with an 
excitations for 

of the FW for 
to be 



M 
Y 



g0Q^12 
4.^UJlUJ2 ' 



(32) 



which is refered to as the theoretical gain in the main 
text. 



For experiments on FW it is important to know also 
the width of the instability regions. The boundaries 
of the unstable region can be found from inspection 
of Eq. (31), which shows that at the boundary we have 
from N that hi = 0^12 /(4<^i) and correspondingly from 
M we then have ai = i^\/^i/^2 5 such that the fre- 
quences to linear order in £0 obtained from (26) are 



2 I ^12 



2 I ^12 /^2 
^^02±^0^W^ 



(33) 

Analogously the system (14) for the case of driven 



boson-boson interactions can be investigated and the re- 
sults are reported in Eqs. (21), (22) and (23) in the main 



text. Below we sketch the derivation also for this case. 

For £0 = the upper equation in (14) is a Mathieu 
equation 



Vu + H 

with solutions V = 
a = 4cJi/(7^ and q 
tation of the cosine- 
We now look again 



■abujl cos{Qt)]V = 0, 



(34) 



A Ce(a,(7,t) + B Se(a,g',t), where 
= —2ai)UJi/Q'^ in the standard no- 
and sine Mathieu functions j34] . 
for solutions to (14) of the form 



of Eqs. (26) and (27). In the case of 0^5 <C 1 one can 
use the expansions Ce{a^q^t) ~ cos(a;oit) + 0{q) and 
Se(a, g,t) ~ sin(a;oi^) + ^{q) Eq. (|34|. In particu- 
lar we then have that the system in the upper equations 
of (28), and hence (29), applies also here and we have in 



the linear order 



n,tt ■ 



UJ02Q1 



-2(5o,tT - ^iQo ^SoVq. 



(35) 



Hence, Aq and Bq are not coupled to Co and I^o in the 
lowest order in eq^ such that it is enough to consider the 
two envelope functions Aq and Bq in the remainder. We 
then require the absence of the resonant terms in the 



right hand side of the upper equation in ( 35 ) . With the 
ansatz Q = 2uji and UJ2 = oJi {bi = ai), we obtain a 
system of equations for Aq (with Aq^t ^ pAo/so) and Bq 
such that the characteristic equation give 



go 
2uji 



Hence the theoretical gain for ai is 



4cji 



(36) 



(37) 



Since p becomes imaginary when > 0^5/4, we have 
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